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NOMENCLATURE 


Symbol Description 



e 

Mt 

p 

0 i 

Subscripts 

c 

P 

o 


Superscripts 
( )’ 

C) 

( ) 


acceleration of gravity 

particle concentration 

drag coefficient 

parameters in turbulence model 

injector diameter 

drop diameter 

mixture fraction 

particle mass flux 

turbulence kinetic energy 

dissipation length scale 

particle mass 

number of particle groups 

number of particles per unit time in group i 

radial distance 

Reynolds number 

source term 

particle source term 

time 

eddy lifetime 

axial velocity 

particle velocity vector 

radial velocity 

Favre radial velocity 

axial distance 

particle position vector 

lateral distance of particle disp^. • o 

path length of particles in an eddy 

time of particle residence in an eddy 

rate of dissipation of turbulence kinetic energy 

turbulent viscosity 

density 

turbulent Prandtl/Schmidt number 
generic property 

centerline quantity 
particle property 
injector exit condition 
ambient condition 

fluctuating quantity 
time mean value 
vector quantity 


ii 
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SUMMARY 


A theoretical and experimental study of particle-laden jets 
is described. The objective of the work is to assist the 
development of 3pray models considering multiphase flows which 
are readily made monodisperse and avoid problems of drop 
shattering and coalescence. Models of these processes, developed 
during the first phase of this study [1-5], were evaluated by 
comparison of predictions with the measurements completed during 
this investigation, a3 well as measurements of Elghobashi et al. 
Lo f 7J. Three models of the process were evaluated: (1) a 

locally homogeneous flow (LHF) model, where slip between the 
phases was neglected; (2) a deterministic separated flow (DSF) 
model, where slip was considered but effects of particle 
dispersion by the turbulence were ignored; and (3) a stochastic 
separated flow (SSF) model, where effects of interphase slip and 
turbulent dispersion were considered using random sampling for 
turbulence properties in conjunction with random— walk 
computations for particle motion. All three models used a k-e 
turbulence model which was extensively evaluated for constant and 
variable density single-phase jets during earlier work in this 
laboratory. 


Measurements were conducted in particle-laden jets, to 
supplement existing results in the literature. Particle Sauter 
mean diameters of 79, 119 and 207 pm, and loading ratios of 0.2 
and 0.66 were considered. Measurements included mean and 
fluctuating velocities of both phases, particle mass fluxes, 
particle size distributions, and calibration of particle drag 
properties. Particular attention was given to defining initial 
conditions of the flows, since the absence of this information 
was a major limitation in obtaining definite model evaluation 
when using existing results. 


The LHF and DSF models did not provide very satisfactory 
predictions for the present measurements. The DSF model 
generally underestimated the rate of spread of the dispersed 

r!? aS fup 3 3 result of ignoring effects of turbulent dispersion. 

The LHF model provided reasonably good predictions for flows 
containing tracer-like particles, but was unsatisfactory for most 
practical flows. The LHF model generally overestimated the rate 
of spread of dispersed phases due to neglect of slip. The SSF 
model, however, yielded satisfactory predictions of flow 
structure— except at high particle loadings. 


Some effects of particles on turbulence properties 
(turbulence modulation) were observed at high particle mass 
loadings. A modified k-e turbulence model, accounting for the 
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effects of interphase transport, was developed for particle-laden 
jets within the framework of SSF formulation. Computations using 
this model showed that the effects of turbulence modulation on 
the present measurments were small — on the order of 7 % at most. 
This finding was expected, since the present test flows are 
relatively dilute. 

A sensitivity study was conducted for predictions of the 
present measurements, to examine the effects of specification of 
initial conditions and particle properties on the calculated 
results. In general, the predictions of gas-phase properties 
were not sensitive to specifications of initial conditions, while 
variations of particle size exerted large effects on predictions 
of particle flow properties. Since the particle size 
distributions involved in this study had relatively small 
standard deviations, uncertainties due to the use of one average 
particle size, the Sauter mean diameter, for each size group were 
within experimental accuracy. 

The particle- laden Jet structure measurements by Elghobashi 
et al. [6,7] were also examined during the present investigation. 
It was found that axial pressure gradients and local 
recirculation zones probably were present in the flows studied by 
these authors. Since the information concerning the magnitude of 
pressure gradients and the extent of recirculation zones are not 
available, these measurements provided only for a qualitative 
evaluation of the models. Nevertheless, the models provided 
reasonable predictions of these measurements. 


1 . INTRODUCTION 


Vi/ 




The objective of this investigation was to complete 
measurements of the structure of particle-laden jets, useful for 
evaluation of models of the process. The experiments considered 
dilute solid-particle-laden jets in a still environment. This 
arrangement has simple geometry and well-defined boundary 
conditions, which facilitates model evaluation. Furthermore, 
dilute solid-particle-laden jets highlight effects of particle 
dispersion by turbulence in jets, while minimizing complications 
due to particle collisions, interphase heat and mass transfer, 
effects of adjacent particles on interphase transport rates and 
modifications of turbulence properties by particles. 

During the first phase of this study [ 1 - 5 ]* models of 
nonevaporating sprays and particle-laden jets were developed and 
evaluated, using the existing measurements in particle-laden 
jets, as well as the new measurements in nonevaporating sprays. 
While predictions for some models were encouraging, the 
evaluation was inhibited throughout by insufficient information 
concerning initial conditions. Moreover, the flow structure data 
available in the existing measurements of particle-laden Jets 
were not complete in most of the cases, which further hindered 
the model evaluation. 

The present investigation supplements the data base of 
existing measurements, considering particle-laden jets involving 
three particle-size groups and two loading ratios. The injector 
properties were carefully characterized so that the measurements 
can be employed to evaluate separated flow models. The flow 
structure measurements included: mean and fluctuating velocities 

of both gas and solid phases, gas-phase Reynolds stress, particle 
mass fluxes and the calibration of particle drag. 

During the course of this study, new measurements in 
particle-laden jets were reported by Elghobashi and his coworkers 
[o,7j. In contrast to all the other existing measurements, these 
authors provided some information on initial conditions of their 
test flows. Profiles of mean and fluctuating velocities of both 
phases and the Reynolds stress of gas phase were reported at x/d 
20. Limited information concerning particle concentration 
profiles was also reported. 


The present new measurements, as well as measurements of 
Elghobashi et al. [6,7], were employed to evaluate particle-laden 
jet models developed earlier in this study. Since these flows 
only involve interphase momentum exchange, modeling efforts 

^Numbers in brackets denote references. 
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concentrated on effects of turbulent fluctuations on momentum 
transfer between the phases, as well as the dispersion of the 
particulate phase by turbulent fluctuations. 

The structure of sprays and other particle-laden flows is 
generally influenced by turbulent dispersion of the discrete 
phase. Turbulent dispersion of particles is examined during this 
investigation by comparing predictions of several theoretical 
models with existing measurements. The study is limited to 
solid-particle-laden jets in a still environment. These results 
are also of interest for spray modeling, however, since the 
geometry approximates near-injector conditions while 
consideration of solid particles avoids complications due to 
polydisperse drop size distributions and drop coalescence. 

Past models of turbulent particle-laden jets often consider 
two limiting cases instead of treating turbulent particle 
dispersion L 5 j. At one limit the particles and the continuous 
phase are assumed to have equal rates of turbulent diffusion. 

The locally homogeneous flow (LHF) approximation provides a 
consistent formulation of this limit. This implies that 
interphase transport rates are infinitely fast, so that both 
phases have the same velocity at each point in the flow. The LHF 
approximation provides best results for flows containing small 
particles, where characteristic response times of particles are 
small in comparison to characteristic times of turbulent 
fluctuations. LHF models have been extensively evaluated during 
earlier work in this laboratory, but only yielded accurate 
predictions for particle sizes smaller than most practical 
applications [ 1 -5 , 8 - 10 ] . 


Turbulent particle dispersion is neglected entirely at the 
other limit. This implies that particles follow deterministic 
trajectories since they only interact with mean properties of the 
continuous phase, yielding a deterministic separated flow (DSF) 
model. Such an approximation is appropriate for flows containing 
large particles, where characteristic particle response times to 
1 low disturbances are large in comparison to characteristic 
turbulent fluctuation times. Several spray and 
solid-particle-laden flow models have been proposed along these 
iines, e.g., El Banhawy and Whitelaw [ll], Mongia and Smith [12], 
and Boyson and Swithenbank [13], and Faeth and coworkers [1-5I 
among others [5]. DSF models generally underestimate the rate’ of 
flow development and particle dispersion for both sprays and 
particle-laden jets, due to the neglect of turbulent dispersion. 


Most practical particle-laden flows exhibit properties 
between these limits and require consideration of turbulent 
particle dispersion. Early dispersion models, discussed by Yuu 
et al. 1 . 1 4 J and Dukowicz [l 5 ] , apply a gradient diffusion 
approximation with empirical correlations for turbulent particle 
exchange coefficients. This approach is not practical, however 
since such exchange coefficients are influenced by both particle 
and turbulence properties requiring excessive effort to 
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accumulate a data base sufficient for general application of the 


Several recent studies of turbulent particle dispersion use 
stochastic separated flow (SSF) methods to circumvent the 
limitations of the gradient diffusion approach. Stochastic 
analysis requires an estimate of the mean and turbulent 
properties of the continuous phase. Particle trajectories are 
then computed by random sampling to find instantaneous continuous 
phase properties. Mean and fluctuating particle properties are 
found by Monte Carlo methods — where a statistically significant 
number of particle trajectories are averaged to obtain system 
properties. 

SSF models have been applied to particle-laden jets. Yuu et 
al. 1.1 h J use empirical correlations of mean and turbulent 
properties for SSF analysis of their particle dispersion 
measurements. Gosman and Ioannides f 1 6 ] propose a more 
comprehensive approach, where flow properties for the stochastic 
calculations are computed with a k-e turbulence model. This 
approach has been adopted by the present authors in their study 
of particle-laden flows, after only minor modification [1-5]. 

In the following, the models are described first of all 
followed by the discussion of experimental methoas used for’ the 
measurements in particle-laden jets. The report concludes with 
an evaluation of the models using the present new measurements 
and measurements of Elghobashi et al. [6,7], along with the 
discussion of difficulties encountered in the experiments of 
Elghobashi et al. [6,7]. 


2 . THEORY 

2.1 General Description 


The theoretical model considers 


o a steady, axi symmetric, 
dilute solid particle-laden, turbulent jet in an infinite, 

stagnant media. It is assumed that the boundary layer* 
approximations apply for the present flow. A k-e turbulence 
model is used to find the properties of the continuous phase, 
3 ince thi 3 approach has been thoroughly calibrated for both 
constant and variable density jets [8-10,17] and has moderate 
computation requirements. The injector exit Mach number is less 
than 0.3; therefore, kinetic energy and viscous dissipation of 

33 Wel1 as 8 as_de nsity variations, are neglected 
with little error. 


2.2 Locally Homogeneous Flow Model 


The formulation of the LHF model corresponds to the 
treatm f nt of the continuous phase for all three models. 
The LHF approximation implies that both phases have the same 
ve.ocity at each point in the flow, i.e., local thermodynamic 
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(2.4) 
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Eq. (2.4). In this instance, it is not necessary to solve a 
transport equation for mixture fraction fluctuations and to adopt 
a probability density function (PDF) for f — which is usually 
necessary during LHF analysis when heat and mass transfer effects 
are considered [5,8-10]. 

Ambient values of u, f, k and e are all zero for the flows 
treated here. Gradients of these quantities are also zero at the 
axis. 

It was generally necessary to approximate initial conditions 
for existing measurements. When slug flow is assumed at the 
injector exit, all properties were taken to be constant, except 
for a shear layer having a thickness equal to 1 percent of the 
injector radius at the passage wall. The constant property 
portion of the flow was specified as follows: 

x - 0 , r < 0.99d/2; 



f - 1 , k - (0.02u ) 2 

o o o 


(2.5) 


e - 2.84 x 10" 5 u 3 /d 
0 o 

Equation (2.5) provides the inner boundary conditions of the 
shear layer until it reaches the jet axis. The initial variation 
of u and f is taken to be linear in the shear layer. Initial 
values of k and e in the shear layer were found by solving their 
transport equations while neglecting convection and diffusion 
terms. 


When fully-developed pipe flow wa£ assumed at the injector 
exit, f 0 was taken to be unity while u 0 was obtained from power 
law expressions provided by Schlicting [21 ] — allowing for 
variations with Reynolds number. Initial values of k Q and e 0 
were obtained from Hinze [22] for fully- developed pipe flow in 
the present Reynolds number range. 

2.3 Deterministic Separated Flow Model 

Both separated flow models adopt the main features of 
the LHF model for the continuous phase. The density of the 
continuous phase is constant, for the present flow, simplifying 
Eq. (2.1). Furthermore, a solution for f is no longer needed, 
since particle concentration is found from the discrete-phase 
solution. 


The dispersed phase is treated by solving Lagrangian 
equations of motion for the particles. At the initial condition, 
the particles are divided into n groups (defined by initial 
position, size and velocity) to yield a statistically significant 
representation of dispersed-pha3e properties. 
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where n, is the number of particles per unit t me in each group, 
"in” and "out” denote conditions entering and leaving a 
computational cell, and Vj is the volume of computational cell J. 

Source terms should also appear 1 " ths Boverning actions 
for k and e-representing direct contributions of the relative 
mnYi nn nf the oarticles to the production and dissipation of 

turbulence. These terms represent the tur V^! nce F ^ d ^i^e n 
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The disperse-phase properties are gained by ~neral form of 

^Basset, Boussines, and Oseen), after neglecting effects of 
particle rotation, can be written as follows 124 J . 
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2 * ** Stochastic Separated Flow Model 

„ taMaM SS f m ° del involves finding trajectories of a 
statistically significant sample of individual particles as thev 

InJe<,t0r “ d e " countei - a ranL dUtriStJr 

of turbulent eddies— using Monte Carlo methods. The treatment of 
the continuous phase is identical to the DSF model. 

Key ele " lents in the SSF model are the method used to specify 
eddy properties and the time of interaction of a particle and a 

fon™^ 7 edd> ’ The approac £ u ? ed to flnd theseVoperties 

follows Gosman and Ioannides [l6], but differs in some details 
roperties are assumed to be uniform within each eddy and to 
change randomly from one eddy to the next. Particle trajectory 

Eqs PU (2 9rLid r ( e 2 t 1 1 0i Sa h e ** ^ ° SF mode1 ' lnv »lving solution of 
tqs. 2.9) and (2.10); however, mean-gas properties in these 

equations are replaced by the instantaneous properties o? each 

nArtJU? - r HS ePti f S ° f ® aCh 6<ldy are f0Und at the Start Of 

orobahnif dd H in ff radtion *>y making a random selection from the 
probability density function (PDF) of velocity. Velocity 

fluctuations are assumed to be isotropic with a Gaussian PDF 

v hav ng a standard deviation of (2k/3)l /2 and mean components u, 

’ °* T ^ e cumulatlve distribution function for each velocity 

llTctZ ^n t T ted , and 3a " pled - ™ 3 

8 three nurnbers in the range 0-1 and computing the 

di«?frih y M 0niP £ nent3 f ° r eaCh three values of the cumulative 
^i^ butl °J faction. This procedure provides a random 

?lictuaUons. " Pr ° pertleS which satisfies the PDF for velocity 

. . t par ^ cle . is assumed to interact with an eddy for a time 
which is the minimum of either the eddy lifetime or the transit 
time required for the particle to cross the eddy. These Umes 

the Till ? a3 ® uraing that the characteristic size of an eddy is 
the dissipation length scale, which is given as [l6,2?] 7 


L e = k 3/2 /c 

Eddy lifetime is computed using 
t e - L e /(2k/3) 1/2 


( 2 . 12 ) 


(2.13) 


a :f assumed to in: act with an eddy as long as the 

. . interaction, At, and the distance of interaction I Ax I 
satisfy the following criteria interaction |Ax p |, 


At < t 
— e 


:™s; 

the interaction la endec with the aecond criterion. 
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The remaining computations are similar to the DSF model, 
except that more particle trajectories must be considered to 
obtain statistically significant particle properties (generally 
2000 trajectories were used). A by-product of the additional 
calculations, however, is that the SSF model yields both mean and 
fluctuating particle properties. This provides an additional 
test of model performance. 

2. 5 Turbulence Modulation Model 

The models discussed in the preceding sections ignore 
direct contributions of interphase transport on turbulence 
properties, i.e., turbulence production, dissipation and scale. 
The comparison between predictions and measurements suggests that 
this approach is acceptable in dilute, lightly- loaded flows. 
However, effects potentially due to turbulence modulation can be 
significant when particle loading ratios are large [5-7,28-31 ], 
requiring modification of the turbulence model to account for the 
direct contributions to turbulence properties due to interactions 
with the discrete phase. 

Several two-phase versions of turbulence models have been 
reported; however, they are limited in application mainly due to 
the assumption of a linear drag law [ 32 - 35 ]. The introduction of 
several new empirical constants which must be calibrated is 
another difficulty with most of the existing turbulence 
modulation models. 

In contrast, the SSF formulation provides a convenient means 
for treating the nonlinear inter phase transport with minimal 
empiricism. Therefore, a modified k-e turbulence model 
accounting for effects of turbulence modulation is derived in the 
following for particle-laden jets within the framework of the SSF 
formulation. 

Since only free jets in stagnant air are considered in this 
study, the pressure gradient is neglected throughout the 
derivation. Following the procedures described by Hinze [ 22 ], 
Tennekes and Lumley [ 36 ], and Daily and Harlow [37], the exact k 
and e equations are obtained from the instantaneous continuous 
phase momentum equations — considering the particle source terras. 
The exact equations are then approximated by the boundary layer 
assumptions. Closure of these time-averaged equations is 
achieved by modeling the turbulent correlations, following the 
work of Gosman, Lockwood and Syed [l 8 ] . The conventional k-e 
model equations are recovered with only one added term due to the 
contribution of the discrete phase in each of the modeled k and e 
equations. The detailed derivation can be found in Appendix A. 
The final form of the equations is as follows 
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l C p Ge)+ i^ (r p v -o e , .I^ (r ^|£, 


ter^ £o7h c °"J aln ® d ln th » b roken-llne boxes" l-l t'he added 
erms due to the interphase transport, and S DU is the particle 

source term in the streamwise momentum equatfon for the 
continuous phase. 

formuiatioi nS ?^ ta HH°H S / 10W properties are known in the SSF 
fpZ ’ 6 term in the k -equation is exact. The new 

term in the e equation, however, is modeled in a way consistent 

tZT i0 T k : £ "° del3 - ° nly or ' e " a “ 00„s?ani SL 

calibration ’ WmCh greatly reduces the effort of model 

3- experimental methods 
3.1 Introduction 

. . . . y he experimental apparatus was used for measurements in 

during Jh^oresenf 3 "? parti cle-laden jets. The flows examined 
during the present study are as follows: 

1 . An isothermal air jet. 

2. Four particle-laden jets (with Sauter mean diameters 79 
119 and 207 um and loading ratios 0.2 and 0.66). 

A variety of measurements were performed to define the fiou 
Ax al and radial profiles of mean and fluctuating ve?ocn?es 
both phases, gas-phase Reynolds stress, and mean 8 particL^f " 
flux were measured. Initial conditions were obtained at x/d - 1 

P article a drag m properties^ eC * PartlCle 3126 di3tributip " a a " d 




3.2 Test Apparatus 


13 


'£>• 


The requirements for a monodisperse particle-laden gas 
jet having no zones of recirculation, a simple geometry and 
well-defined boundary conditions can be met with the apparatus 
illustrated in Fig. 1. The injection tube was mounted on a 
two-dimensional traversing mechanism at the center of the cage. 
For all the test flows, injection was downward into 3 till room 
air. The test cage was 1 m square by 2.5 m high— enclosed with 
1 b -mesh screens to reduce the influence of room disturbances. 
Traversing in the third dimension, to obtain radial profiles of 
flow quantities, involved moving the entire cage assembly, which 
was mounted on a bearing track. This configuration allows a 
fixed optical measurement station, minimizing alignment problems 
for the laser-Doppler anemometer ( LDA) system. 

Measurements from this arrangement should be attractive for 
those wishing to evaluate models. The flow is parabolic; 
therefore, its turbulence characteristics can be modeled more 
reliably than flows with recirculation. Also, parabolic flow 
greatly simplifies problems of obtaining accurate numerical 
solutions. Boundary conditions are well-defined, since there are 
no uncertainties regarding wall friction and inlet flow 
properties which are encountered in confined chambers. 

The flow system is illustrated in Figure 2. The Jet tube has 
an internal diameter of 10.9 mm and extends in the vertical 
direction for 90 injector diameters. Flow at the exit of the 
injector roughly corresponds to fully-developed turbulent pipe 

flow; however, the initial condition i 3 completely measured in 
any event. 

The air supply for the particle-laden jet apparatus was 
provided by an oil-free air compressor. The air flow was metered 
using a critical flow orifice. Seeding particles needed for LDA 
measurements of ga3-phase quantities were added using a 
reverse-cyclone seeder. Sand particles for two-phase flows were 
introduced upstream of the injection pipe by a gear feeder 
powered by a variable speed motor. The particle mean flow rates 
were varied by changing motor speed, as determined by system 
calibration. Three size ranges of particles, with Sauter mean 
diameters 79, 119 and 207 pm, were used in the experiments. The 
sand was sifted to reduce the diameter variation within each size 
group. Typical size distributions are illustrated in Figure 3 
where more than 1000 particles were sampled for each group. The 
density of the sand particles is 2620 kg/m3. 


3« 3 Experimental Procedure 

3*3.1 Gas-Phase Velocity Measurements 

Table 2 is a summary of the instrumentation used 
i or each measurements. 
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Particle size distributions for the particle-laden jet experiments 







't)* 


C 

U J3 H 

B 

u * 
a) o M 
0,10 co 
co ^ h 


u 

a> cd 

03 44 

OJ >> CO 
H U-O 
, G 
m a) ju 
<0 3 (0 
•u cr 4J 

*2 ?> g 
tfl H 3 
O M-i o 
co a 

#» 

T3 ^ 4-1 
M CO CO 
nj o ^ 

> g 3 

M 0 

O O 

** S4= 

a) 4J • 
* G ’H GO 
g <0 S G 
CO tH 

0 M T3 co 

JQ d) <u CO 

1 H 4J QJ 
H 0,4-4 O 
CO CX *H O 

3 o x: u 

Q Q co q, 


a) 

cc 

to i u 
X c o 
a, o 44 
i a 

co *4 

co a <u 

GO *rl 44 

* 3 

- a, 
owe 

4-4 4-JO 

q o 

CO <0 *H 
CO Q PI 

a) *h • 

CO M 0 00 

0 3 c 

3 co jr *h 
cr co 4 -j co 
•rt <o tI co 
C B £ (D 
X o 

o >> q o 

dJ 44 O Vr 


<0 0 O CO 
6 »—♦ q 44 
(0 <D 3 cd 
w > nTJ 


rH S p O 

O C C 

•H d) co 


4-1 CO 
CO 

HD -rl 

a -u 

CO -H CO 
a co 
G O d> 
cfl H M 
<0 <0 44 

se > to 


GO 

G CO 
•H <0 
4-J -rt 
CO 44 

9 

rH 

44 

G 

(0 

•H 

U 

•rt 

4h 

44 

(0 

o 

3 *H 

44 

u 

•u a 
o o 

CO 

GO 

3 iH 

CO 

cd 

» — 1 <D 

CO 

U 

4-i > 

6 

O 

T3 CD 

CD 

CO 

G rH 

rH 

rH 

CO U 

O 

O 


*H 

*H 

G 44 

44 

44 

CO >4 

Wl 


CD CO 

CO 

3 

as a. 

Pu, 

P-, 


Mean and fluctuating gas velocities were measured using a 
s mgl e-charnel, dual-beam, frequency-shifted LDA. An equipment 
list for the LDA system appears in Table 3. The sending and 
receiving optics had a focal length of 2 42 mm with a 11.61° angle 
between beams. The signal was focused on the photomultiplier 

W * 3 i?° ? m focal len Sth lens. The aperture diameter of the 
photomultiplier was 0.25 mm. The receiving lens was masked with 
a beam-stop which provided a collection aperture of 25 mm 
diameter. The preceding optical configuration produced an 
ellipsoidal measuring volume of 0.470 mm in length and 0.098 mm 
in diameter, with a fringe spacing of 3.128 urn. 

The injected flow was seeded with 0.2 um A1 2 0o particles, 
which was sufficient for most of the measurements? This primary 
seeding was supplemented by seeding the surroundings when 
measurements were taken near the edge of the jets, where the 
primary seeding level was low, in order to avoid biasing. 

axial^d^f^ 0 ^ 11 ! 1 ! veloclty °°®P°nents were measured in the 
axial and radial directions at several axial distances from the 

injector for all test conditions. The use of several beam 

°ri!" fcati0nS allowed measurement of Reynolds stress. Reynolds 
stress measurements are valuable for checking model predictions 

f T ?° Catlng the flow centerline (since Reynolds stress is 
zero at the centerline and doesn’t exhibit the broad maxima 
encountered for other variables) . Concentration biasing and 
interference from large particles were avoided by employing high 
concentrations of seeding particles. An integration time of 
one minute or more was used to determine the time-averaged 

quantities. Each measurement was repeated three times to get a 
better averaging. 8 


3*3.2 Particle Velocity Measurements 

, n « P artlcle ve lodty was measured using the same 

eain a sp?M ement r aS f ° r gas ~ phaae measurements. In this case, the 
T °f * P hoton,ulti Pli^ and counter were adjusted to 
elatively low levels so that only strong scattering signals from 
large particles were recorded. The output of the data processor' 

oolleoted " lth a 11/23 minicomputer and processed 
to yield mean and fluctuating particle velocities (by taking 
number-averaged velocities). When these tests were run, no 

seeding was used for the gas phase, in order to reduce potential 
bias errors of the measurements. potential 

po „ h Th ® ^ enfc f rline Particle velocity for all the test flows was 
echecked using the double-flash imaging technique developed by 
a. 5. p. Solomon in a companion study. The photograph system 

focus 3 the ?L£r° 3 f" mlcro3econd flash sources, a lens system to 
focus the light and a camera. The flashes were fired 

consecutively at electronically controlled times so that two 

images of the moving particle, were obtained cn the aa^e 

photograph. Measurement of the distance travelled by the 
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Table 3 LDA Equipment List 


Component 

Manufacturer 

Model 

Helium-Neon Laser 

Spectra Physics 

125 A 

Integrated Optics 

Thermo-Systems 

900 

Frequency Shifter 

Thermo- Systems 

9180 

Photomultiplier 

The rmo- Sy s t ems 

960 

Counter 

The rmo-Sys t ems 

1980 

RMS Voltmeter 

Thermo- Systems 

1060 

Dual Beam Oscilloscope 

Tektronix 

561 A 

Integrating Digital 
Voltmeter 

Hewlett-Packard 

240 IC 



at' 


TllllloT^lTocnT. lnteriral 6etuee " th * t» flashes provided 


3 ' 3 * 3 Particle Flux Measurements 

sampling probe, ^he^iobe^a^a^odinprt^^ ?* th a " i30kineti c 

collection probe used by Szekely and FaethTlfil 0 " ri ^ particle 
probe had an ID of 3 mm uhi^h Y . , Faeth L 38 J - The tip of the 

resolution for the me^^ents^^ ^f Sufficlent 3 P a tial 
maintained by applying suction tr> ^° kine J ic conditions were 
sampling flow rate. The samDlimr ri^ Pr<>be 30(1 meterin « the 

ZotT tM to — 

a balance which^r^sSsitim^^ Wei8hed wlth 

times varied from 1 minut* ti on f ? ne miligra, “- Sampling 
location in the^low ^! narM n ^ e3 ” depending on «» 
the cross-section of the jet at each » U f data ” inte «rated over 
against the injector particle mas* ri * lal position_ “ w as checked 
sampling efficiency? P3rticle mass flow r *te to determine the 


3 . 3 . 


Measurements 


the shape of the^amcles^as solewhaf T** necessar y since 
particles might not follow thf °® e ” hat Irregular and the 
spheres. These raeaqurwnn * ? standard drag law for solid 

velocities of the particles^n^ree-fair 6 ^ 1 " 1 ?? terninal 
gas velocities resulting from cartel f a’ 1 33 the local 

drag coefficients were computed usin* U3lng th * LDA ’ The 

the Sauter mean diameter of each si z f IrouT^ Vel ° Cities 30(1 

were'obtained^by^adj ustin^the^standard Heyn0lds nurabers 

of measured drag to that of standard dr f g law » usir,g the ratio 
Reynolds number for each Darti^? ™ at the free - fa H 
numbers Involved in the dra» m» a group * The Particle Reynolds 
from i to 10, c^pieS to tho;! 6 ” 6 ^ 3 Were sma11 ’ ranging 
jet measurements, where the Dart ®”?° uatered in the particle-laden 
the range of 1 ti 200 depend?^ nUmbera variad in 

the particle size. However the sam» 6 po f) tion in the flow and 
over the heynolds nunher ra^ge USed 

** • RESULTS and discussion 
4 * 1 Experimental Cnn/nn^ 


oondUtoJ^rT;^^^ ^enTs^' I single-phase 63 6 


(Pr 
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a i r jet and four dilute particle-laden jets, 
particle size groups and two loading ratios, 
f u test conditions are summarized in 

test results are tabulated in Appendix B. 


involving three 
were chosen for 
Table 4 and all 


foundry ?and ln - the experiments “ere Penn State C-30 
roundry sand, which was sieved to obtain different size ranges 

k^/m 3 a The n^t 31 ^ ° f the particles was measured to be 2620 ' 
k ? h' ^ Th P 3 ^ 1016 dra 8 properties were calibrated against th* 
standard drag law for spheres. A multiplication factor Jo the 
standard drag law was found for each of "the site groJp" u! 

207’ rn’ 5 and !< 51 fOP partlcle groups having SMD of 79 , 119 a nc 

IZ The fact that the wild", do not" bey 

tne drag law for spheres exactly is reasonable due to the 
somewhat irregular shape of the partlcles--partloularl5 for 

Xr- The oaii&rated *■« ia -= — irs/* 

let I«nt ^? i mea3ureraent3 of mean Petiole velocities along the 
Jerhnf 1 were checked U3in 8 the double-flash photographic 

reasonable flT * greement between the two techniques was 

° WS containin 8 sma H and intermediate size 

case Ttet<, e ’ 8 u' WUhin 5% f ° r 0336 1 jet and 15? for ca 3® 2 and 
were sicmif particle velocities measured using LDj 
signif icantly lower than those measured usins double-f lA^h 
Photographic technique for case 4 particle-laden !e?7sm - S? 

x "; f° 2 Tto r 1?»°at X/ /d = ^ -gld'fr™ 'to ? a, 

n A. Tki 1 V* 1 /d * 3 °’ where the LDA measurements yielded 
Lot Li 3Up The exact cause for Hi “ 

PDA measurements for large particles Is not clear. Due t^thla 

»»locitiee h f LD# ” ea f urements ° f "«an and fluctuating particle 
oaae 4 particle-laden jet are not employed f^r 
-valuation of predictions in the following. 


^•2 Initial Conditions 


1 for all fi vl fi if 3 f b0th P hases w ere measured at x/d 

Thy, iat . 1 f V f 1 ows . The results are presented in Figures 4-1 
The jets expanded well beyond r/x = 0 B th* 8 ^? 3 H ' 

lot IT !*' T P ° nent or veloclt y fluctuation, „as 

S2 ^'enie^^^c^er;^ 31 *° ' al °“ ia “‘"< 

.f 11 bhe gas-phase results, except near the jet edge were 
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Summary of Test Conditions 
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Fig. 7. The initial conditions for case 2 particle 

laden jet radial profiles of gas-phase n» 

and turbulent quantities at x/d = 1. 
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Fig. 8. The initial conditions for case 2 particle-laden 
J et radial profiles of solid-phase mean and 
turbulent quantities at x/d = 1. 
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Fig, 9. The initial conditions for case 3 particle 

laden jet radial profiles of gas-phase 

mean and turbulent quantities at x/d * 1- 















The presence of particles resulted in no significant change 
for the ga 3 -phase properties at x/d * 1. The mean velocity 
profiles were very similar for all the five test flows. However, 
the turbulence kinetic energies in the particle-laden jets were 
slightly higher than in the single-phase jet. 

The profiles of mean particle mass flux, mean and fluctuating 
particle velocities and particle turbulence kinetic energy are 
shown in Figures 6, 8, 10 and 12. The mean slip velocity ratios, 
(u-u D )/u, at the jet axis were 7.4$, 18.9$, 13.2$ and 26.8$ for 
case 1 to case 4 par tide- laden jets, respectively. The slip 
decreased with loading ratio for the jets having the same 
particle size (compare the case 2 and case 3 flows) which agrees 
with measurements in Refs. 6 and 7. The mean particle velocity 
profiles were wider than the mean-gas velocities, except for the 
case 4 flow due to the difficulties discussed previously. 

The particle mass fluxes were nearly uniform up to r/x ■ 0.5 
for particle-laden jets with larger particles, and then decreased 
to zero at edge, cf.. Figures 8, 10 and 12. The jet with 
smallest particles, however, exhibited a continuous variation of 
particle-mass flux over the injector exit, cf.. Figure 6. 
Particle- turbulent quantities were slightly smaller than 
comparable gas-phase quantities--except for the case 4 flow as 
discussed earlier. 



4.3 Axial Variation of Flow Properties ‘ 

Predictions of the SSF model will not include effects of 
turbulence modulation for results illustrated in Sections 4.3 end 
4.4, since these effects are not large for the present test flows 
and we wish to consider model performance with all empirical 
factors established earlier. Turbulence modulation will be 
considered separately in Section 4.5. 

Predicted and measured profiles of centerline mean gas-phase 
velocity along the axis are illustrated in Figure 13. 

Predictions for the single-phase jet, based on the model of 
Shearer et al. [8] developed earlier in this laboratory, agree 
well with measurements. This establishes an acceptable baseline 
for measurements in the particle-laden jets, since the 
s ingle- phase jet model was well-calibrated [5,8,9j. The measured 

rates of decay of centerline velocity in the particle-laden jets j 

>■ are smaller than in the single-phase jet, with this effect more 

l! pronounced at higher loading ratio. Predictions of both the LHF | 

f and SSF models for particle-laden jets are in fair agreement with ' 

| measurements, except for case 3 flow where both theories , 

overestimate the rate of centerline velocity decay. j 

] 

Predicted and measured mean-particle velocities along the Jet j 

axis are illustrated in Fig. 14. The SSF model predictions are j 

in good agreement with the measurements in all cases. The LHF j 

model, however, significantly overestimates the rate of velocity j 

l - 





0/ 


OF POOR QO. .LiTV 



Fl8 ' ° f CM “ rll " C ««" solid- phase 


35 



V 


decay at all conditions, since the test particles have 
significant inertia. 

The axial variation of gas-phase turbulence-kinetic energy is 
illustrated in Fig. 15 . Results are shown for both 
particle-laden jets and the air jet. The theories agree well 
with measurements except in the near-injector region, where all 
the models underestimate turbulence levels. The underestimation 
of turbulence kinetic energy near the injector exit could be due 
to limitation of the k-e model, since the turbulence is still 
developing in this region. The fact that turbulence levels in 
fully-developed regions, approaching x/d = 40, roughly correspond 
to values estimated from models which ignore effects of particle 
motion on turbulence properties, indicates that effects of 
turbulence modulation were small in these flows. This is 
reasonable, since all the test flows were dilute. 

The axial variation of centerline particle turbulence kinetic 
energy is illustrated in Fig. 16. The results for the case 4 
particle- laden jet are not shown, due to problems with the LDA 
measurements of large parti cl es—discussed previously. The LHF 
model significantly overestimates the particle turbulence levels 
for all the cases considered here, since the particles cannot 
actually follow the gas-phase turbulent motion due to their 
inertia. In contrast, the SSF model yields fair agreement with 
measurements. 

Predicted and measured particle-mass fluxe 3 along the jet 
axis for the four particle-laden jets are illustrated in Fig. 17. 
The LHF model significantly overestimates the particle 
dispersion, due to neglect of effects of slip between the phases. 
In contrast, the SSF model provides satisfactory predictions in 
all cases. 


4. 4 Radial Variation of Flow Properties 

Theoretical predictions of the radial variation of mean 
and fluctuating quantities for both phases were compared with 
measurements completed during the present study. The comparison 
involved profiles of mean velocity, axial and radial components 
of fluctuating velocity, turbulence kinetic energy and Reynolds 
stress for the gas phase; and profiles of mean velocity and mass 
flux, axial and radial components of fluctuating velocity, and 
turbulence kinetic energy for the particle phase. The tangential 
component of flu ctua ti ng ve locity was not measured; therefore, it 
was assumed that v'^ = w^ — which is usually true for Jets and 
sprays [28]— to find measured turbulence kinetic energies for 
both phases. The predictions of gas-phas e fluctua ting velocity 
components were obtained by assuming u’2 ;v i 2 „ 1:0>5 k--which 
is usually observed in single-phase jets. 
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Predictions and measurements for air jet are compared in 
Figs. 18 and 19. Predictions agree well with measurements for 
all quantities, just as in the case of Shearer [8] and Mao [9]. 

Predictions and measurements for case 1 particle-laden jet 
are compared in Figs. 20-23* In this case, the loading ratio was 
quite low; therefore, gas-phase profiles approach the air-jet 
properties. Both LHF and SSF predictions are in good agreement 
with measurements for gas-phase properties at x/d = 20, except 
that the SSF model underestimates Reynolds stress at its maximum. 
The agreement between predicted and measured ga3-phase properties 
is less satisfactory at x/d = 40. The predicted radial profiles 
are generally too wide compared to experimental results at this 
axial position. Predictions of particle properties are compared 
with measurements in Figs. 21 and 23. SSF predictions are in 
good agreement with data for all properties, except for the axial 
component of fluctuating velocity, where the model predictions 
are too low near the axis and too high away from the axis. 
Predictions of the LHF model are not satisfactory in this case. 
Results for the DSF model are also shown in Fig. 23. The DSF 
model yields poor results, similar to the findings when 
evaluating these models using existing measurements. Neglecting 
particle dispersion by turbulence causes the rate of spread of 
the flow to be substantially underestimated. 

Comparison of radial profiles of mean and fluctuating 
quantities for the case 2 particle-laden jet is shown in Figs. 
24-27. In this case, effects of particles on gas-phase 
properties are small due to the low loading ratio and the larger 
particle size. Predictions of gas-phase properties for both the 
LHF and SSF models are in reasonably good agreement with 
measurements at x/d = 20 and 40. The particle properties are 
also well predicted by the SSF model. In contrast, the LHF model 
overestimates all particle properties while the DSF model 
underestimates the flow spreading rate — similar to the case 1 
flow. 

The results for case 3 particle-laden jet are illustrated in 
Figs. 28-31. The loading ratio in this case is higher, resulting 
in greater effects of particles on the properties of the gas 
phase. The measured mean-gas velocity profiles for the case 3 
jet are narrower than the case 2 jet, which has the same particle 
size, but a lower loading ratio. The fluctuating velocities, as 
well as the Reynolds stress, also decrease as the loading ratio 
increases. This may be attributed to the effects of turbulence 
modulation, i.e., turbulence intensity levels are attenuated by 
particle drag, resulting in a reduced jet spreading rate. This 
trend cannot be reproduced by the models, since effects of 
particle drag on turbulence structure are ignored for these 
predictions. As noted earlier, effects of turbulence modulation 
on the present test flows will be discussed in Section 4.5. 

The particle properties are well predicted by SSF model, 
except for the profile of particle mass flux at x/d « 40 (Fig. 
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Fig. 30. Radial variation of gas-phase mean and 

turbulent quantities in case 3 particle- 
laden jet at x/d = 40. 
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31 ) where the SSF model somewhat overestimates the 
parti cle-spreading rate. The fact that good agreement is 
obtained for particle properties, in spite of errors in 
predictions for gas phase, indicate the particles in this test 
fLow are relatively insensitive to the variation of gas 

properties due to their large inertia (large size and initial 
velocity) • 

Predictions of gas-phase properties for the case 4 
particle-laden jet are compared with measurements in Figs. 32 and 
33- The results for particle properties are not shown, since in 
this case velocity measurements for large particles using the LDA 
experienced difficulties. At this test condition, the large 
particles have a relatively small velocity change and short 
residence time from the jet exit to the measuring stations. 

Hence the influence of particles on gas-phase structure is small, 
even though the loading ratio is high. The predictions of the 
SSF model are in good agreement with measurements for all 
gas-phase properties, while the LHF model slightly overestimates 
both mean and fluctuating quantities. 

In general, the SSF model provided reasonably good 
predictions of mean and turbulent quantities for both solid and 
gas phases. This is encouraging, in view of the simple 
formulation for treating the interphase transport processes. 
Although the LHF model yields poor predictions for solid-phase 
properties, it is capable of providing reasonable estimation for 
gas-phase properties of the flows examined in this study. 

However, it should be recalled that the present measurements were 
limited to dilute flows with relatively large particles, where 
the effects of interphase transport on gas-phase structure are 
not very significant. The DSF model is not attractive for the 
present flows, in view of the poor performance in predicting 
particle properties and the need for detailed information for 
initial conditions. 



5 Effects of Turbulence Modulation 

For results thus far, the SSF model has not considered 
direct contributions of interphase transport on turbulence 
properties. The comparison between predictions and measurements 
reveals that this approach is acceptable in dilute flows. 
However, effects of turbulence modulation can be large when 
particle-loading ratios were large, e.g., the results of Laats 
and Frishman [ij and the present study in Figs. 28 and 30 . 

Effects of turbulence modulation were studied for the 
measurements completed in the present investigation, using the 


'•'-O'* f wo X 1 IK 

model discussed in Section 2.5. The new model constant, Co 
which should be of order unity based on physical reasoning, * was 
allowed to vary between 1 and 5. The calculated results for the 
fractional change of centerline flow properties, due to effects 

of turbulence modulation, are presented in Table 5 , for Co 

o*1 
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Fig. 33. Radial variation of gas-phase mean and 
turbulent quantities in case 4 particle- 
laden jet at x/d = 40. 



Summary of Results of the Turbulence Modulation Study 





























7K 

and 5. The maximum change of flow properties from this 
calculation is 6 . 6 ?, which occurs in turbulence kinetic energy at 
i x/d = 20 for case 3 particle-laden jet. Furthermore, the 

' predictions were found to be relatively insensitive to the 

I variations of C 3 . 

| The small effect of particles on the turbulence structure of 

j presently measured flows is expected, since these flows are all 

relatively dilute. The large particle inertia and small particle 
residence time may al 3 o reduce the effects of interphase 
transport on turbulence properties for these flows. 

The relatively large discrepancies between predictions and 
measurements, shown in Figs. 28 and 30 cannot be rectified by 
I only considering effects of turbulence modulation in the 

predictions. This deviation between theory and experiment is 
more likely to be due to experimental errors. 



In the preceding discussion, it was noted that 
specification of initial conditions is of vital importance to 
predictions using separated flow models. Because of this 
concern, special care was exercised to obtain reliable 
measurements of initial conditions for the test flows during the 
present experimental investigation. However, uncertainties still 
exist in the specification of particle properties, i.e., the 
particle-drag coefficient and mean-particle size used for 
initiating the computations. As discussed previously, the 
particle-drag coefficient for each size group used in the 
predictions was obtained by adjusting the standard drag law with 
a single scaling factor obtained from the free-fall experiment. 
Some uncertainties were introduced in this process, since the 
scaling factor for each particle group might vary over the 
Reynolds number range experienced in par tide- laden jet 
experiments. 


The particle size is another concern. The measured 
particle-size distributions for the three size groups have 
standard deviations varying from 11$, for largest particle-size 
group, to 23$, for smallest particle-size group. Only one 
average size, i.e., the Sauter mean diameter, was used for each 
group— to save computational effort. 


Another potential source of uncertainty in predictions comes 
from the specification of e Q , the rate of dissipation of 
turbulence kinetic energy at the jet exit. The k-e model 
requires initial values of k and e to start the computation. 

WhUe k 0 was measured directly. e 0 was calculated from measured 
values of turbulence kinetic energy, Reynolds stress and the mean 
velocity gradient at exit, using Equation (2.3) and the gradient 
tansport assumption. Although this approach was consistent with 
the formulation of k-e model, errors might be introduced due to 
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the small magnitude of the velocity gradient near the jet 
centerline at the exit. 

Due to these uncertainties regarding specification of initial 
conditions and particle properties, the sensitivity of model 
predictions to the variation of these parameters was examined. 

The results of this study are presented in Tables 6, 7, 8 and 9 
for the four particle-laden jets investigated in the present 
experiments. The entries show the fractional change in 
predictions brought about by raising the value of the indicated 
parameter by 1 00% — with all other parameters unchanged. These 
estimates were obtained using the SSF model without considering 
effects of turbulence modulation. The entries are computed for 
the major measurements stations considered in the 
experiments — x/d - 20 and 40. 

In general, the predicted flow properties for both gas and 
particulate phases, are not sensitive to the variation of e 0 and 
k 0 . This is expected, since the production and dissipation of 
turbulence kinetic energy are the prevailing mechanisms for the 
development of turbulence structure in free Jets. In contrast, 
the calculated solid-phase properties show larger sensitivities 
to the variations of particle diameter and drag coefficient. 

This finding also illustrates the importance of defining initial 
particle properties in order to obtain a convincing evaluation of 
models. 

Broadly, the most influential parameter for ttie predictions 
is particle size. For example, a 51% change in Up C and -50% 
change in kp/u^Q was obtained by raising the value of d p by 100%. 
The large sensitivity of d p is due to the fact that particle 
relaxation time (the time for the particle to change the slip 
velocity by a factor of e“l) is roughly proportional to d p ^. 

Since the standard deviations of the particle-size distributions 
involved in the present study are small (varying from 11% to 23 % 
for the three size groups), the errors introduced by assuming 
monodisperse particle- laden flow are within experimental 
accuracy. 


5. MODEL EVALUATION: MEASUREMENTS OF ELGHOBASHI ET AL. 


5. 1 General Description 

Recently, Elghobashi and his coworkers [ 6 , 7 ] reported 
LDA measurements for particle-laden jets, with particle sizes of 
50 pm and 200 pm and loading ratios from 0.32 to 0.85. The most 
important contribution of their experiments is that, in contrast 
to all the other existing measurements, they provided some 
information on initial conditions of their test flows. The 
authors also paid special attention to eliminating the 
cross-interference between the two phases in the velocity 
measurements. Flow structure measurements included: axial decay 

of streamwise mean velocities along the jet axis for both phases. 
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Table 6 

Summary of Results of Sensitivity Study for 
Case 1 Particle-Laden Jet a 


f 

U 





Output Variables 

— 

Input 

Parameter 

u 

c 

k /u 2 
c c 

u 

pc 

_ 2 
k /u 
pc pc 

G 

c 

x/d = 20 

e 

o 

.09 

-.01 

.07 

-.12 

.06 

k 

0 

-.12 

.08 

-.06 

.09 

-.06 

C D 

.05 

.02 

-.13 

.22 

-.12 

d 

P 

.04 

-.03 

.18 

-.40 

.13 

x/d = 40 

1 , 

0 

.05 

'.00 

.09 

-.11 

.08 

k 

0 

-.06 

.02 

-.08 

.13 

-.08 

C D 

.03 

.04 

-.30 

.28 

-.16 

d 

P 

.06 

-.04 

.51 

-.50 

.18 


Entries show the fractional change in 
about by raising the indicated input 


prediction brought 
parameter by 100%. 


» 



Table 7 

Summary of Results of Sensitivity Study for Case 
2 Particle-Laden Jet a 


Input 

Parameter 


x/d = 20 


Output Variables 


k/u 2 
c c 


k /Z 2 

pc pc 


£ 

O 

.09 

-.02 

.02 

-.03 

.06 

k 

o 

-.13 

.05 

-.06 

.04 

-.04 

C D 

.02 

-.01 

-.09 

.14 

-.09 

d 

P 

.02 

-.02 

.06 

-.21 

.08 

x/d = 40 






£ 

O 

.04 

~.00 

.04 

-.04 

.06 

k 

0 

-.07 

.01 

-.09 

.07 

-.07 

C D 

.02 

.01 

-.32 

.19 

-.12 

d 

P 

.02 

-.03 

.27 

-.26 

.13 


Sh °? ! he f 5‘ aCtlonal chan 8 e In prediction brought 
about by raising the indicated input parameter by 100%. 
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Table 8 

Summary of Results of Sensitivity Study for Case 
3 Particle-Laden Jet a 



Input 

Parameter 


Output Variables 


u 

c 

k /u ^ u k /u ^ 

c c pc pc pc 

G 

c 


x/d = 20 






£ 

o 

.08 

-.02 

.02 

-.04 

.06 

k 

o 

-.13 

.04 

-.07 

.05 

-.05 

C D 

.04 

-.02 

-.10 

.12 

-.09 

d 

P 

x/d = 40 

.03 

-.02 

.05 

-.18 

.08 

£ 

O 

.03 

-.00 

.06 

-.06 

.07 

k 

o 

-.06 

.01 

-.05 

.06 

-.07 

C D 

.05 

.03 

-.31 

.19 

-.14 

d 

P 

.05 

-.02 

.30 

-.26 

.16 


3 

Entries show the fractional change in prediction brought 
about by raising the indicated input parameter by 100%. 



Table 9 


Summary of Results of Sensitivity Study for Case 
4 Particle-Laden Jet a 


Output Variables 


Input 


Parameter 

u 

c 

k/u 
c c 

_x/d » 20 

£ 

O 

.09 

-.02 

k 

o 

-.14 

.04 

C D 

.02 

-.02 

d 

P 

.01 

-.02 

x/d » 40 

£ 

O 

.04 

-.00 

k 

o 

-.07 

'.00 

C D 

.04 

.02 

d 

p 

.02 

-.02 


u 

pc 

k/u 2 
pc pc 

G 

c 

.01 

-.02 

.03 

1 

o 

.03 

-.03 

1 

• 

o 

O' 

.09 

-.06 

.01 

-.10 

.05 


.02 

-.03 

.03 

-.03 

.05 

-.04 

-.22 

.15 

-.08 

.18 

-.14 

.07 


Entries show the fractional change in 
about by raising the indicated input 


prediction brought 
parameter by 100%. 
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adial profiles of mean and fluctuating velocities for both 
phases and Reynolds stress of the gas phase, at x/d = 20. 

Limited information concerning particle concentration profiles 
was also reported. 

^nJ^ rae f S f^f tS ° f Elghobashi et al. are unusually complete 
1 Cat J d f LM and data Processing systems were used to obtain 

mwn f 1 ondlUons at the J'et exit, as well as flow properties at 
downstream positions. However careful examination of the test 
apparatus employed in their experiment indicated that axial 
pressure gradients and local recirculation zones were very likely 
ThP^ pr ® 3ent the flows investigated by Elghobashi et al. 

a^d altered'the f/^T 6113 AOt defined in the measurements 
and aitered the flow structure and complicated the specification 

of boundary conditions for predictions. 

of a T ri e ^ COnfl r atIOn 0f E1 8^l=a3hi, et al. [6,7] consisted 
a 2 om Jet issuing into a cylindrical gla 33 chamber with yerv 
leu velocity coflou. The chamber diameter uas 60 cm. The jet 

u^re'o 0°? mL a Ld X ^,“^ 'l'* the <*>«« veSclues 

,, . * 3 and 0,1 m/s for test cases, respectively. All 

located It the^otj 3 ^ 8 ^ through the exhausfc ducting system 
located at the bottom of the test chamber. This arrangement 

feedi^ ex f reme f care for operation of the exhaust and coflow 
eding systems to balance the inflow and outflow to the test 
chamber, and to match the entrainment flow for the jet 

florfllld^whn 1 f l0W ,L ieldS an 3Xial P ressure gradient over the 
field, while insufficient supply of coflow air generates 
recirculation zones near the jet boundary. 

5 « 2 Results and Discussions 

The experiment of Elghobashi et al. [6,7] provided 
measurements in single-phase air jets, in addition to The 
measurements in particle-laden jets. This air jet data was fir<*i 

thiq l ? t h l U f lnS the single_ P hase J efc model developed earlier in 

fJe oa« f8 ZY’ Mi" 06 th ‘ 3 "° del ^ been «t«slvely tested In 
p l >9J. No recirculation zone was assumed in the 

boundary condition specifications for predictions. The measured 
and predicted centeriine velocity decay are summarized in ?ab^e 

that withon^L^ 3 U6S ° f 3X131 pressure gradient. It is evident 
fiPTH ^ 1 ad ? ing a . constant pressure gradient in the flow 
model significantly overestimates the rate of 
centerlme velocity decay. It is noted that the same model 
yielded satisfactory results over varieties of test conditions 
for^flows similar to Elghobashi et al. [6,7] in the past 

exhanif C n!! ai ? tain J nR balance betw een total inlet flow and 
exnaust flow is extremely difficult, it is believed that axial 

pressure gradients were present in the test chamber. In order to 
in more insight a study was conducted to investigate the 
potential presence of recirculation zones in the How n e Jd it 

is known that wnen the coflow falls short of the entrainment 
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capacity. of the jet in a confined duct, the balance is satisfied 
wi recirculated fluid- The entrainment calculation using the 
cmpm.cal expression of Abramovich [no], showed that all the 
coflow air was entrained at x/d - 1 4 and x/d - 26 for the let. 

This f?nir alr ' ,el30ltles of °-°5 m/s and 0.1 m/s, respectively 
This finding suggested that recirculation zones were present In 
he regions where measurements were conducted. The results of 
the entrainment calculation was checked by testing for the 

Beo"e e r Ce H: t:. r a e „ C rSnf U ° n ??,"? ' JSln8 tha t«hnl q 8 ue described by 
rhP Wiliiams [41 . The calculated front edge of 

e s. 

extend 

From the preceding discussion, it is evident due tn fh<» 
arrangement of test apparatus, the flow field Investigated bv 
Elghobashl et al. [6,7] was subjected to the InflSnofVVaVL 
pressure gradients and local recirculation zones. sLetS 

Vnd°tw atl °V °° nc * rnln 6 the magnitude of the pressure gradients 
and the extent of recirculation zones were not reported bv the 

M/mf inal authors * two Pressure gradients, i.e., 0 N/m3 and - 5 0 
’ were c ° n3 idered in the predictions discussed in the 

InfoS™. reclr ® la “°" “ere Ignored, for lack of 

h h t, Pr ® dlCted and measur ed radial profiles of flow proDerties fnr 
both phases at x/d - 20 are Illustrated In Fig. 3 J ?pr^te 

ratio oTo a v" "T a 3 ''!? 8 3 paPtIcU of 50 um and a loading 

ignoring Jhe iiffeVts^ the SSF “ ,oael - both and 

SSF ? e [ fects of turbulence modulation, denoted SSFM and 

SSF respectively, are shown on the plot. The mean velocity 

^n^i e V re " onnali2ed wit * the centerline velocity of the 
° ngle phase jet, to illustrate the magnitude of the sli D 
velocity and the effects of particles on the gas ptase Velocity 

is #> clear frora the figure that improvement in the Y 
predictions of the SSF model is obtained by adding a -5 0 N/m3 
pressure gradient, which further suggests Ihe presence ^ a 

gradient in the flow field. The SSFM model generally 
yields better results than the SSF model for most of ?he now 

e?Ve P «s U of VVib 3 f rated ln F1S - nis fact Indi^atVs th^ 

effects of turbulence modulation are important in this test flow 
The overall agreement between SSFM model predictions and 
measurements is fair, except for axial pamc!rflu«uaUon 
velocity, where the model undereatlmates the fluctuation level. 

It is interesting to note that the effects of narUnniskm 

lareer° n ^ ^ turbulent Properties of thecas phase are 

larger in this case than those shown in Fig. 20 for the 

s 0 m“nVr b »r? a T° le ’ 1, - :en JCt - 13 <>« to the fonowlngi 

lnjwur em°£o fifl partIcle fasldence time from ?he 

larger-- half f he ™ easurln « Nation (roughly four times 
g alf the initial particle velocity and twice the axial 
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distance) , and the potential effects of the axial pressure 
gradient and recirculation zone, for this test condition. 

Predicted and measured radial profiles of flow properties for 
both phases at x/d = 20 are compared in Fig. 35 for the 
particle-laden jet with 50 urn particles and 0.85 loading ratio. 

I he agreement between model predictions and measurements is 
simitar to the results with the lower loading ratio, illustrated 
in Fig. 34. However, as particle loading increases, the rate of 
decay of the gas-phase centerline velocity is reduced. This 
trend is reproduced by the SSFM model. 

Results fc- larger particles, 200 pm in diameter, and loading 
ratio of 0.64 are illustrated in Fig. 36. Agreement between 
measurements and predictions from the SSFM and SSF models, with 
-5.0 N/md pressure gradient, is satisfactory for all the flow 
properties, except for the mean particle velocity profile, where 
both models underestimate the particle velocity far from the jet 
axis. In this flow, effects of turbulence modulation are not 
important, since the large-size particles are less effective for 
interphase transport. 

Agreement between predictions and measurements is generally 
better for the 200 pm particle jet than for the 50 ym particle 
► ul-,.* maln difference between these flows is that, besides 
p i^ ference ln P article 3ize » the coflow velocity is 0.1 m/s 
for 200 ym particle jet and 0.05 m/s for 50 ym particle jets. 

This is interesting, since the recirculation zone, from the 
preceding discussion, extends to upstream of measuring station 
ix/d 20) for the case with coflow velocity of 0.05 m/s. The 
front boundary of recirculation zone only reaches to slightly 
downstream of the measuring station for the case with 0.1 m/s 
coflow velocity. Therefore, the less satisfactory agreement for 
the cases with 50 ym particles may be due to effects of the 
recirculation flow penetrating upstream of the measuring station 
since this phenomenon is not considered in the present model. 

6- SUMMARY AND CONCLUSIONS 
6. 1 Summary 

The overall objective of the present study was to 
develop and evaluate models of dilute turbulent particle-laden 
jets. Three models of these flows were developed: (1) a 

locally homogeneous flow (LHF) model, where interphase transport 
rates are assumed to be infinitely fast; (2) a deterministic 
eparated flow (DSF) model, which allows for finite interphase 
ransport rates (evaluated using mean properties of the 
continuous phase) , but ignores dispersion of the particle Dhase 
by turbulent fluctuations; and (3) a stochastic s'eparaied^S 
SSF) model, where finite interphase transport rates, and the 
turbulent dispersion of particles are considered by allowing 
particles to interact with a succession of eddies using a 
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random-walk procedure. Typical models of each type were 
described and evaluated using measurements in dilute 
particle-laden duct flows and round jets. All versions of these 
models used a well-calibrated k-e model to estimate mean and 
fluctuating properties of the continuous phase. 

The theoretical description of the continuous phase for all 
three models was based upon the Reynolds-averaged form of the 
conservation equations written in Eulerian coordinates. The 
dispersed phase is treated, for the separated flow models, by 
solving Lagrangian equations of motion for the particles. A 
modified version of the GENMIX program [42] combined with a 
second order Runge-Kutta ordinary differential equation solver 
for particle motion were used to solve the governing equations. 

Experiments were conducted for a single-phase air jet and 

fh^ P Mf rti f le lad n n ^? ts ’ to su PPlement existing measurements in 
the literature. Particle Sauter mean diameters of 79, 119 and 
207 um, and loading ratios of 0.2 and 0.66 were considered. 
Measurements included mean and fluctuating velocities of both 
phases, parti cie mass fluxes, particle-size distributions, and 
calibration of particle-drag properties. A single-channel, 

DA WaS ^ t0 obtain velocity and drag measurements, 
articles were examined with a microscope to find size 

distributions. The particle-laden flows were isokinetically 
sampled to obtain profiles of particle-mass flux. These 
measurements were compared with predictions of all three models 
discussed in this study. 

A modified k-e model, incorporating direct contributions of 
interphase transport on turbulence properties (turbulence 
mo u ation), was developed within the framework of the SSF model. 
The potential effects of turbulence modulation on the present 
measurements were examined using this model. 

A sensitivity study was conducted to investigate the 
in luence of uncertainties in specification of initial conditions 
and parti cie properties of model predictions. The results are 
useful in identifying potential sources of error for both 
predict! ons an d meas ur emen t s . 

™? a3Urement3 . COmpleted by E1 Shobashi , et al. present 
ifficulties in obtaining a definitive evaluation of models. The 

3b ™“ d w®* apo ^ nd the experiment was the source of most of the 
difficulties. In the first place, the shroud causes an axial 
pressure gradient of unknown magnitude. Secondly, the coflow of 
?n i^r 3 n- n3U f fiCl ^ nt t0 prevent recirculation of the flow beyond 
, 7 ™ the Jet exlt ‘ Based on comparison between 

we nd tha 3 t a ^ mea ? U, : en,ent3 f0r their 5l ^le-phase experiments 
nf o!r d M / h f t ^ ial Pressure gradient should be on the order 
N/m * The two-phase flow calculations were completed 
enng pressure gradients in the range of 0 to -5 N/m3 but 
ignoring the effect of the recirculation pattern. Calculations 
were also conducted both including and ignoring the effects o? 




turbulence modulation. In view of the uncertainty in 
experimental conditions, the agreement between predictions and 
measurements is reasonably good. 

6. 2 Conclusions 

The major conclusions and observations of this study are 
as follows: J 


T. The LHF and DSF models did not provide very satisfactory 
predictions for the data base considered in this 
evaluation. The LHF model was only satisfactory for 
flows containing tracer-like particles, where 
characteristic response times of the particles are small 
in comparison to all characteristic response times of the 
continuous phase. The DSF model generally underestimated 
the rate of flow development and particle dispersion over 
the entire data base. Both of these methods appear to 
have limited utility for modeling practical 
particle-laden flows. 

2. The SSF model yielded reasonably good results over the 
entire data base (both existing and present 
measurements)— with no modification in the prescription 
for eddy properties from its original calibration. This 
result is encouraging, since the SSF method has the 
potential to treat nonlinear and complex interactions 
between the phases which are encountered in practical 
particle-laden flows, e.g., evaporating and combusting 
sprays, condensing and reacting vapor jets in liquid, 
pulverized coal combustion, etc. Additional development 
and evaluation of the SSF method will be required to 
determine if this potential can be realized. 

3. Effects of turbulence modulation were small for the flows 
investigated in the present experiments, since these 
particle-laden flows were all lightly-loaded with 
particles having relatively large 3 izes and short 
residence times-which tends to inhibit interphase 
transport. However, results presented in Table 5 
indicate the effects of turbulence modulation can be 
significant for applications involving higher loading 
ratios and smaller particle sizes than the present test 

l UWo * 

4. In general, present model predictions are relatively 
insensitive to the specification of gas-phase initial 
conditions. The specification of particle properties 
however, exerts much more pronounced effects on 
predictions. This emphasizes the importance of 
measurements of initial par tide- phase conditions and 

^ cle pr °P ert ^ es » drag and size, in order to 

obtain a convincing evaluation of separated flow models. 
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applied to a wide ran^P nr n f' these models can be 
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appendix a 


DERIVATIO N OF k-EQUATION AND 

WITH PRESENCE OF DISPERSED Phacp 


A. 1 Derivation of k-Eguatlon 

Since free jets issuing into st 
i-tudy, the pressure gradient is nec 


agnant air are considered in this 


' SraGlent 1S neglected throughout the following 

derivation. The governing equations are limited to steady, axisvn^etric, 
constant property flows. 


The gas-phase momentum 


equation, including the parti cl 


e source term. 


) =Vr — A-+ F 
j 1 3 dx j° X j P u j 


(A. 1) 


where F pu ^ - s^/o, end the repeated indices imply summation. 
Multiplying E,. (A.l) by „ t and rearranging terms yields 


~ 1Uj + 

J J 1 T A 


J J 


F (A. 2) 
pu i ' 


T " e 1 ”“ n “ n ““ «« ^ decomposed to mean and fluctuating 

components according to the following expressions 


U. * u. + u.’ 
1 i i 


Vi = “i “i + 2 u ± u ± ' + q : 


s k 


(A. 3) 
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After decomposing and Reynolds-averaging, Eq. (A. 2) may t 


written as 


"sT <1 Vi» 


o , (1 2 \ 6 - rl Z'] dr- 

' “j t? 1 ) - 5JJ “j )• SJ K “i u j 


12i 


+v 


Dx 


Du 3 u 9 u . . ~ 2 Du! Du! 

- (u — . v _i - v I ia_ , V _JL i 

j 1 ^ x j Dxj dXj 2 Dxj 5xj Dx^ 


+ u. F 
l pu. 


(A. 4) 


The Reynolds-averaged form of Eq. (A. 2) is 


_D 

9x 


3u . 


: ( Vp ■ 57T <V 5T - “1 “p + f p„ w-5) 

J J J i 


Multiplying Eq. (A. 5) by and rearranging terms yields 


D - A - - 


dU, 


3u. 3u. 


.. - — u.(-ru.u.)] = V T (u 7 — — ) - V r — — - 

JX .s J 2 ii 3x. i <ix. Ox, Ox . 

J J J J J 


5u i 3 


+ u i u j 5T - 5T ( “i u i “P + “i F 


P“< 


(A. 6) 


Subtracting Eq. (A. 6) from Eq. (A. A) the equation for turbulent kinet 
energy is obtained 
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.» ,2 


J J J 1 J 3x, 3x.3x. 


t'u| Du! 

v ( sir ^ (aT^) + u. F 


"x 3x. 

J J 


■i r - u. F 
1 P u • i pu . 


(A. 7) 


The viscous diffusion term in Ea f A 7\ • , 

q» (A. 7) is much smaller than 

U '” 0t " er temS ^ be ne8lecte ' i - The fourth ter. in the RHS 

cue equation the viscous dissipation of turbulent kinetic 


energy, c . 


After employing the Deling approach of Gasman et al. ( 18 ), the 
modeled thin shaar layer for. of *. (A . 7) C a„ be wrlt[e „ ^ 

v & 2 . c 


+ u F - u f 

Pu pu 


(A. 8) 


where the (v F - Z v \ „ , , — - 

pv py ) nd (w F pw ) terras have been neglected, 
since they are small compared to (G~F -G p ) 

Pu pu ^ * 

Equation (A. 8, is derived for constant density f lous £or 

simplicity and to highlight the derivation of the new terms due to 
interphase transport. however, the derivative of turbulent kinetic 

difficulties. 
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Following past work [18,39], the following equation is obtained 

for variable density flow, after completion of the order of magnitude 
estimation 


5 

dX 


(p u k) + i f- (r p v° k) = i i- (r 
r dr r dr 


|jt 9k 

o, 9r ; 
k 


+ U (0) - P £ + u S 

r or pu 


- u S 


(A. 9) 


where p v - P v + pV and S pu - p F pu . In arriving at this equation, 
a number of density fluctuation terms are ignored and other terms are 
modeled according to the methods employed for uniform property flows. 
In fact, the only major difference between Eq. (A. 9) and the constant 
density form, Eq. (A. 8), involves retention of the p’v' correlation 
term on the LHS of the equation. The justification for this largely 
rests on the success of the method during earlier study of variable 
density flows [8,9,17,18,39]. 


A. 2 Derivation of e-Equation 

The e-equation, considering effects of interphase transport, was 

derived in the same manner as the k-equation. However, since the 

algebra is quite involved, only the procedure is outlined here. The 

first step in the derivation is to obtain a transport equation for the 

fluctuating velocity, u’, by subtracting Eq. (A. 5) from Zq. (A. 1). 

The e-equation is obtained by differentiating the u! equation with 

9 U : 1 

respect to x R , multiplying throughout by v — - and finally time 

D 

averaging. The order of magnitude analysis is then conducted and 


i 

■i 

I 

| 

f 

1 

| 

1 

£ 

4 

J 
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higher order terms are neglected, following the procedure by Tennekes 
and Lumley [36 ] . The resulting equation becomes 


& 


_ ( 1 

3 u . 3u ! du\ 3u, T 3u, ? -u’ ju* ju! 

- (I. E) - - 2v r-i — i r-J- + — - 2V -i r-L 

i dx. dx. ijx, ex. , x. jx, ex, ex. 

J J t «- k J -• xj kkj 


a 2 , .2 , 

d u. d u. 


3u! SF! 


r— u ! C ’ - 2V ' i i — + 2v ^ (A. 10) 

” x j J ax j 3x k Jx j ox k JX k 3x k 


Equation (A. 10) is then transformed to thin shear layer form in 
Cartesian coordinates and modeled according to Gosman, Lockwood and 
Syed [ 18] . The modeled e-equation can be written as 


1 3 / - . 13 t 3c e liu 

“ (u l) + - — (r V e) = - T— (r — — ) + v - — 


^"1 ^ ^ -J t | 


r 2 „ , oF* „ , dF' ,, , 3F' 

_ c L_ + o v £u + EZ EZ + 3w pw 

e k “ dr 3r 3r dr 'dr dr 


(A. 11) 


The last term in the RHS of the above equation is the contribution 
from interphase transport and needs to be modeled. Employing the same 
method for single-phase flow [12], this terra is modeled as follows 


, , 3F 1 . , 3P' a , 3F* 

9 3u' pu + 3v' pv 3w pw 

3r 3r dr dr dr 3r 


" 2v C c 3 ^ < A ‘ 12 > 


where C is the new model constant to be determined. 
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T>,e variable density version of modeled ^-equation is obtained 
following the same argument as for k-equation, as follows 


5 - kh (r^|i 


. £ ,' ju . 2 _ c 2 r :-s 

H t k ( a7> ' c . p T- - 2 c P * 


( a . i n 
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Gas-Phase Properties 
x/d = 1 


r/x 


0.0 

.046 

.138 

.229 

.333 

.413 

.459 

.505 

.551 

.596 

.642 

.688 

.734 


r/x 


0.0 

.023 

.046 

.067 

.092 

.115 

.138 

.161 

.174 


u/u 


r 3 i. 


u’T/u 

fv 2 

i c 1 

^ j 


/u 


1.0 

.986 

.976 

.931 

.888 

.804 

.718 

.621 

.473 

.295 

.215 

.115 

.059 


.063 

.065 

.066 

.067 

.072 

.096 

.120 

.147 

.154 

.136 

.114 

.077 

.058 


.039 

.038 

.036 

.040 

.046 

.058 

.079 

.098 

.114 

.113 

.096 

.075 

.073 


Gas-Phase Properties 
x/d «= 20 


u/u 



C 1 

v’ 2 

» j 

c 

l J 


/u 


1.0 

.927 

.738 

.555 

.392 

.252 

.171 

.119 

.078 


.271 

.270 

.269 

.230 

.197 

.156 

.120 

.084 

.057 


.181 

.182 

.190 

.177 

.156 

.128 

.095 

.070 

.055 
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u’v'/u 


- 2 


.00013 
.00053 
.00052 
.00057 
.00136 
.00226 
.00486 
.00789 
. 01040 
.00887 
. 00602 
. 00384 
.00126 


u'v’/u 


- 2 


.0003 

.0116 

.0159 

.0164 

.0135 

.0098 

.0046 

.0028 

.0008 




E> 













0.0 

1.0 

.113 

.032 

.00017 

.OSS 

.993 

.113 

.033 

.00054 

.156 

.986 

.115 

.038 

.00158 

.248 

.948 

.111 

-043 

-OQUftfr 

.339 

880 

.10 

.054 

.00229 

.431 

.786 

.122, 

.081 

.00365 

.493 

.337 

see ' 
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